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ABSTRACT 


A  new  instrument  for  measuring  the  concentration  and  size  of  raindrops 
at  the  ground  is  described.  This  instrument  measures  the  time  of  flight  of 
drops  across  a  sample  volume  and  thus  does  not  require  a  measured  or  assumed 
drop  velocity  as  do  existing  disdrometers.  The  signals  expected  from  the 
device  are  analysed  and  a  detailed  description  of  the  electronics  complete 
with  circuit  diagrams  is  given.  This  is  followed  by  the  Z8Q  machine  code 
routine  for  measuring  the  time  of  flight  and  size  of  the  droplets  in  real 
time  as  they  cross  the  sample  volume,  together  with  the  high  level  program 
for  subsequently  soring  this  data  and  printing  the  spectra.  The  procedure 
for  setting  up  and  calibrating  the  instrument  is  outlined  and  some  sample 
results  obtained  in  the  field  are  presented.  These  results  show  that  the 
device  is  working  satisfactorily.  Finally  some  suggestions  for  future 
work  are  given. 


Instrument  Summary- 
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The  distrometer  detects  raindrops  as  they  enter  and  leave  a  cylindrical 
sample  volume.  The  sample  volume  is  defined  by  a  sheath  of  light  100  ym 
thick  which  surrounds  the  cylinder.  This  light  is  collected  on  to  a 
photodiode  as  shown  in  Figure  1.  A  drop  entering  and  leaving  the  volume 
will  give  rise  to  two  pulses  from  the  photodiode  amplifier  output  having  height 
proportional  to  the  drop  size  and  separation  in  time  equal  to  the  drop  transit 
time.  As  shown  in  Figure  1  single  pulses  which  result  from  a  drop  passing 
through  the  edge  of  the  annular  beam  are  rejected  from  the  subsequent  computer 
analysis  and  only  pairs  of  pulses  are  accepted.  Thus  for  each  pulse  pair 
we  have  the  amplitude  which  gives  the  diameter  (d)  of  the  drop  and  the  transit 
time  (t). 

a)  Flux  method 

The  first  and  more  conventional  analysis  technique  is  to  calculate  the 
flux  of  drops  and  then  from  an  assumed  velocity  to  find  the  concentration. 

If  the  projected  sample  area  of  the  volume  is  A,  and  in  a  time  T,  there  are 
n  pulses  pairs  corresponding  to  drops  of  diameter  size  d,  then  the  flux  of  these 
drops  is  n/.T.  .  and  if  the  drop  velocity  is  v  then  the  concentration  is 

n/<v.A.T.). 

b)  Time  of  flight  method 

The  seond  method  uses  the  time  of  flight  information.  If  in  a  time  T, 

for  drops  of  size  d,  the  total  time  of  flight  of  all  the  drops  is  It,  we  know 

It  / 

that  drops  of  this  size  were  in  the  sample  volume,  V,  a  fraction  /T  of  the 
time,  and  thus  the  concentration  is  given  as  directly  as  VT/Et  drops  per  unit 
volume.  No  knowledge  or  assumption  of  the  drop  velocity  is  needed,  and 
the  instrument  can  operate  in  conditions  of  high  windspeed. 


2. 
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Comparison  with  existing  commercial  instruments 

The  two  commercial  instruments  available  both  measure  the  flux 
of  raindrops  and  assume  terminal  velocities  to  derive  a  drop  concentration. 

As  shown  below  both  suffer  from  some  shortcomings  in  windy  weather. 

The  new  instrument  described  in  this  report  is  also  of  much  simpler 
construction  and  does  not  suffer  from  any  critical  optical  alignment 
problems . 

a)  Joss-Waldvogel  Distrometer  (1967) 

This  device  measures  the  momentum  of  drops  as  they  impact  upon  a 
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50cm  horizontal  plate.  Earlier  examples  of  this  type  had  an  unacceptable 
dead  time  during  which  the  oscillations  excited  by  each  impact  died  away. 

In  this  device,  however,  a  complex  servo-system  keeps  the  plate  stationary. 

The  minimum  size  detectable  is  generally  300pm  but  this  may  not  be  attained 
during  acoustic  interference  from  strong  wind  or  thunderstorms.  As  the 
sampling  head  is  in  the  form  of  a  toadstooljin  windy  weather  when  the 
drops  are  not  falling  vertically  the  sample  area  will  change  in  an 
unpredictable  manner. 

b)  The  Knollenberg  Optical  Array  Distrometer  (1970,  1972) 

The  PMS  ground-based  optical  array  precipitation  spectra  probe  is 

designed  to  measure  drops  over  the  size  range  0.2  to  12.4mm  with  a  sample 
2 

area  up  to  60cm  .  This  instrument  measures  the  shadows  cast  by  raindrops 
on  a  linear  photodiode  array. 

It  has  been  developed  from  an  aircraft  mounted  instrument  in  which 
case  the  terminal  velocities  are  negligible  compared  with  the  aircraft 
velocity.  The  performance  of  the  instrument  is  described  by  Heymsfield  (1976) 
and  Curry  and  Schemenauer  (1979).  White  Sands  scientists  have  experienced 
difficulties  in  obtaining  reliable  drop  concentrations,  derived  from  the 
drop  flux  measurements  of  the  ground-based  instrument,  during  windy  weather. 


3*  The  UMIST  Disdrometer 


The  present  prototype  device  operates  on  the  shadowgraph  principle,  but, 
whereas  the  PMS  (Knollenberg,  1972)  device  has  a  horizontal  rectangular  sample 
area,  this  instrument  detects  drops  as  they  enter  and  leave  a  cylindrical  volume. 

This  cylindrical  volume  presents  an  equal  area  to  any  drops  which,  due  to  high  winds, 
are  not  falling  vertically. 

Figure  1  shows  the  optical  arrangement.  The  device  is  very  simple  so  that 
operation  in  hazardous  field  environments  will  be  reliable.  The  use  of 
coherent  radiation  has  been  avoided  because  of  interference  problems  arising  from 
dust  particles  and  the  expensive  optical  elements  which  would  be  required.  None 
of  the  lenses  need  to  be  very  precisely  positioned  thus  no  accurate  adjustments 
are  required  in  the  field. 

A  quarz-halogen  filament  light  source  of  diameter  0.5mm  at  the  focus  of  a  lens 
of  focal  length  25cm  provides  a  parallel  beam  of  light.  After  this  light  has 
crossed  the  sample  volume  (in  this  example  a  distance  of  17cm)  it  is  normally 
incident  upon  a  mask  which  only  transmits  light  falling  upon  an  annulus  of  diameter 
3  cm  and  width  lOOym,  thus  defining  the  boundary  of  the  sample  volume.  The  light 
passing  through  the  annulus  is  then  focused  on  a  large  area  photodiode.  When  a 
drop  enters  and  leaves  the  volume  it  should  give  rise  to  two  equally  sized  pulses 
from  the  photodiode  amplifier  as  it  obscures  first  one  side  of  the  annulus  and 
then  the  other.  Provided  that  the  drop  is  larger  than  lOOym  diameter  the  amplitude 
of  the  pulses  should  be  proportional  to  the  drop  diameter.  Single  pulses  resulting 
from  drops  passing  through  the  edge  of  the  annulus  are  rejected  in  subsequent 
computer  analysis.  This  rejection  criterion  is  much  simpler  than  the  linear  array 
of  photodiodes  required  for  the  PMS  instrument, 
a)  Diffraction  limits  and  sampling  volume 

Knollenberg  (1970)  shows  that  a  raindrop  acts  as  an  opaque  sphere  and 
refraction  can  be  neglected  but  the  smallest  detectable  drop  size  is  limited 
by  diffraction.  It  can  be  seen  from  Knollenberg' s  isodensitometer  traces  of 
diffraction  patterns  of  particles  in  white  light  that,  at  a  distance  equal  to 


2  ... 

3d  (where  d  is  diameter  and  X  the  wavelength  of  light)  the  total  intensity 

4X 

of  light  measured  along  a  narrow  slice  through  the  diffraction  pattern  has 
fallen  to  75Z  of  its  value  for  a  perfect  geometric  shadow.  For  a  sampling 
cylinder  of  length  17cm  this  corresponds  to  a  particle  of  diameter  300um  being 
measured  to  an  accuracy  of  lOOum.  Laboratory  calibration  of  the  disdrometer 
with  artificially  produced  drops  confirmed  this  limit. 

2 

The  17cm  sample  length,  giving  a  sample  area  of  51cm  ,  was  chosen  as  a 
compromise  so  that  a  reasonable  number  of  large  drops  would  be  sampled  without 
having  too  high  a  minimum  'detectable  drop  size.  If  it  were  considered  important 
to  build  up  a  statistically  significant  number  of  both  very  small  and  very  large 
drops  in  the  minimum  time  then  two  instruments  with  differing  lengths  of  sampling 
cylinder  could  be  operated  simultaneously.  The  disdrometer  is  constructed  so 
that  this  is  a  trivial  mechanical  adjustment. 


b)  Setting  up  procedure  and  tolerances 

The  apparatus  has  been  designed  to  be  rugged  and  simple,  with  no  precise 
adjustments  required  in  order  to  facilitate  field  operation.  The  0.5mm  diameter 
light  source  placed  at  the  focus  of  a  5cm  diameter,  25cm  focal  length  lens  results 
in  a  beam  divergence  of  1  in  500.  In  order  not  to  exceed  this  divergence, 
due  to  the  finite  source,  the  source  must  be  placed  so  that  it  is  within  2mm 
both  axially  and  radially  of  the  lens'  focus.  Over  17cm  the  1  in  500 
divergence  will  introduce  a  maximum  penumbra  of  340ym,  however,  because  the 
annulus  defines  a  slice  through  the  shadow  and  penumbra  the  error  introduced 
is  much  less  than  this. 

A  large  area  photodiode  (3mm  x  3mm)  is  used  so  that  the  10cm  lens, 
collecting  light  from  the  annulus,  needs  to  be  placed  to  within  3mm  of  the 
axis  and  along  the  axis  to  within  1cm  of  the  focal  point. 

The  only  precision  required  lies  with  the  annulus.  Is  is  most  important 

that  the  width  be  100pm  (+  1  Z)  and  that  the  transmission  around  the  annulus 


be  uniform  to  1Z. 
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c)  Engineering  Drawing 

In  the  appendix  is  a  full  drawing  of  the  optical  arrangement  of  the 
distrometer.  The  light  source  is  a  6V  10W  miniature  quartz-halogen  light 
source  type  M29  manufactured  by  Thorn  Lighting,  U.K.  The  source  is  mounted 
so  that  the  filament  coil  axis  is  along  the  axis  of  the  instrument  to  keep 
the  source  diameter  a  minimum. 

4.  Photodiode  Pulse  Characteristics 

Before  describing  the  electronic  design  and  circuitry  in  detail  it  is 
important  to  calculate  the  various  characteristics  of  the  pulse  expected 
from  the  photodiode  current  to  voltage 
(a)  Individual  pulse  risetimes 

The  smallest  drop  detectable  is  limited  by  diffraction  to  300pm  diameter 
and  thus  a  pulse  length  of  lOysec  would  correspond  to  a  drop  of  the  size 
crossing  the  lOOym  wide  beam  at  30m  sec  Larger  drops  or  lesser  winds 

will  result  in  longer  pulses.  The  longest  pulses  will  be  caused  by  the 
largest  drops  falling  at  terminal  velocity.  For  a  4mm  drop  at  8m  sec  ^ 
the  pulse  length  will  be  500ysecs.  The  electronics  should  be  able  to  cope 
with  pulse  lengths  of  10-500ysecs  or  rise  times  of  half  that  time. 

(b)  Pulse  pair  separation 

The  possible  range  of  separation  of  the  pulse  pairs  is  summarised  in 
Table  1.  The  drops  may  either  fall  at  terminal  velocity  or  at  the  speed  of 
the  wind  (a  maximum  value  of  30ms  ^  being  assumed),  and  they  may  cross  a 
diameter  (30mm)  or  near  the  edge  of  the  beam.  If  the  drop  passes  through 
the  edge  of  the  beam  there  will  of  course  be  only  a  single  pulse.  Assuming 
that  the  two  pulses  will  be  just  resolvable  if  the  light  intensity  returns 
to  normal  between  the  pulses  we  can  use  the  intersecting  chords  theorum  as 
shown  in  Figure  2  to  calculate  the  minimum  resolvable  transit  distance. 
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If  the  cylinder  diameter  is  d  (30mm)  and  the  drop  diameter  is  2x,  then: 

(d-x)x  *  l2 

where  21  is  the  minimum  transit  distance,  and  is  displayed  in  Table  1. 
From  Table  1  for  drops  of  relevant  size  falling  at  terminal  velocity  the 
maximum  transit  time  is  26  msecs  and  the  minimum  is  1.6msecs.  In  the  very 
high  wind  of  30ms  ^  the  pulse  separation  times  range  from  1msec  down  to 
only  140psecs.  If  the  axis  of  the  instrument  is  not  perpendicular  to  the 
wind  then  these  times  will  be  increased  slightly. 


(c)  Probability  of  overlapping  pulse  pairs 

The  probability  of  finding  two  drops  in  the  sample  volume  at  one  time 
can  be  found  from  Poisson  statistics.  If  the  concentration  of  raindrops 
is  n  and  the  sample  volume  is  V,  the  the  average  probability  of  finding 
one  drop  in  the  beam  is  ^  =  n  V.  Poisson  statistics  then  give  the 
probability  of  finding  m  drops  in  the  beam  as : 


It  is  quite  justifiable  to  assume  that  the  drops  are  falling  randomly.  In 
this  case  the  results  are  summarised  in  Table  2  for  a  sample  volume  of 

3 

120cm  (cylinder  diameter  3cm,  length  17cm)  and  assuming  a  Marshall  Palmer 
raindrop  size  distribution.  It  is  well  known  that  this  distribution 
overestimates  the  concentration  for  small  drops  and  we  can  see  that  at  10mm  hr 
the  chance  of  finding  two  drops  greater  than  500ym  in  the  beam  at  one  time 
is  only  one  in  200.  In  the  initial  design  of  the  hardware  and  the  software 
sorting  programs  it  is  proposed  to  ignore  this  probability. 
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5 .  Electronics  and  Circuit  Diagrams 

When  drops  enter  and  leave  the  sample  volume  there  is  a  pulse  on  the 
output  of  the  photodiode  amplifier.  We  need  to  measure  the  size  of  this 
pulse  and  the  elapsed  time  since  the  last  pulse  to  enable  the  amplitude 
of  pulse  pairs  and  transit  times  of  individual  drops  to  be  calculated. 

To  enable  this  pulse  data  to  be  stored  and  then  subsequently  sorted  into 
raindrop  spectra  a  small  Z80  based  'Nascom'  microcomputer  is  used  and  the 
spectra  output  to  a  printer.  The  drops  are  measured  and  timed  using  a 
machine  code  program  in  real  time,  but  the  subsequent  sorting  is  performed 
using  a  slow  program  in  'Basic'. 

The  overall  electronics  block  diagram  is  shown  in  Figure  3  (a). 

The  timing  is  achieved  using  the  Zylog  CTC  chip,  and  as  shown  below,  the 
ports  8,  9,  10  and  11  of  the  Z80  microprocessor  are  connected  to  various 
timers  to  enable  the  drop  transit  time  to  be  found  and  also  the  real  time 
in  minutes  and  seconds  to  be  read  so  that  absolute  concentrations  can  be 
calculated.  The  photodiode  and  pre-amplifier  are  housed  with  the  distrometer 
up  to  20m  from  the  computer.  Finally  we  consider  the  computer  interface, 
which  detects  the  peak  value  of  pulses  and  signals  to  the  computer  via  port  5 
that  it  has  an  eight  bit  pulse  height  ready  on  port  4.  We  shall  now  consider 
these  three  sections  in  detail. 

(a)  The  Counter  Timer  Board 

There  are  two  requirements  for  the  timing  of  the  distrometer.  We  need 
to  measure  the  time  between  the  pulse  pairs  caused  by  a  single  drop,  entering 
and  leaving  the  sample  volume,  which  was  shown  in  Section  4  to  be  in  the  range 
1  to  24  msecs  for  conditions  of  reasonably  low  wind  speed;  and  also  we  need 
to  know  the  total  sample  time  which  will  probably  be  of  the  order  of  minutes. 
Because  of  this  two  external  clocks  are  needed  at  4  khz  and  1  hz. 
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The  CTC  board  has  four  eight* bit  counter  timer  registers,  the  mode  of 

each  counter  timer  being  controlled  by  two  control  words.  The  configuration 

22 

used  is  shown  in  Figure  4.  A  crystal  oscillator  at  2  hz  is  divided  down 
using  4040  counters  to  supply  a  4.096khz  clock  for  the  pulse  separation  timing 
and  a  lhz  for  a  minutes  and  seconds  counter.  The  4khz  clock  is  fed  into 
channel  0  on  the  CTC  board  which  counts  quarter  milliseconds  and  the 
channel  0  overflow  is  fed  into  channel  1.  Thus  channels  0  and  1  (corresponding 

to  Ports  8  and  9)  hold  a  sixteen  bit  time  in  quarter  milliseconds.  The  lhz 
clock  is  fed  to  channel  2,  which  resets  after  sixty  counts,  and  triggers  channel 
3  on  resetting,  so  channels  2  and  3  (Ports  10  and  11),  hold  the  seconds  and 
minutes  counts.  Figure  4  also  shows  the  control  words  needed  to  set  the 
Ports  8,  9,  10  and  11  to  operate  in  this  way. 

(b)  Pre-amplifier  and  line  driver 

The  pre-amplifier  electronics  and  power  supply  are  housed  in  the  same 
diecast  box  as  the  photodiode  and  the  receiving  optics,  and  are  shown  in. 

Figure  5.  The  arrangement  is  capable  of  driving  signals  down  cables  at 
least  20m  long  to  the  computer  with  its  interface  electronics  which  are 
housed  in  some  adjacent  shelter. 

The  overall  circuitry  is  reasonably  straightforward,  but  there  are  some 
practical  points  which  require  discussion.  We  note  that  from  Section  4(a) 
the  pulse  lengths  should  be  in  the  range  10-500ysecs 

The  four  IC's  act  as  voltage  to  current  converter,  unity  gain  buffer, 
a.c.  inverting  amplifier,  and  analogue  line  driver.  Two  outputs  are  available, 
the  output  from  the  buffer  monitors  the  overall  light  level  incident  upon  the 
photodiode  and  is  useful  for  verifying  that  the  alignment  is  correct,  while 
the  pulses  from  the  raindrops  after  further  a.c.  amplification  are  transmitted 
down  the  cable  to  the  computer  interface  using  a  high-current  high-slew-rate 
analogue  driver. 


1 
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The  photodiode  is  a  SD  -  444  -  11  -  11  -  171.  The  50mm  receiving 
area  is  most  useful  in  alleviating  any  alignment  problems  but  the  capacitance 
of  this  area  does  mean  that  any  subsequent  amplifier  is  prone  to  oscillate. 

The  first  stage  current  to  voltage  converter  uses  a  3140  with  1.8M2  and  5pf 
in  the  feedback,  this  gives  a  rise  time  of  9psecs;  any  higher  gain  in  this 
first  stage  leads  to  oscillation,  as  does  a  capacitative  load  for  the  3140. 

As  a  consequence  a  unity  gain  741  buffer  is  inserted  before  the  capacitor  of 
the  a.c.  amplifier  and  the  capacitance  of  the  20m  cable  carrying  the  d.c.  monitor 
level  from  this  current  to  voltage  converter. 

Using  a  light  source  the  6V  10W  miniature  quartz-halogen  source  type  M29 
manufactured  by  Thorn  Lighting  U.K.  mounted  as  described  in  Section  3  the 
monitor  light  level  should  be  about  6V.  The  first  3140  amplifier  saturates 
with  an  output  of  11V  and  consequently  some  slowly  varying  extra  diffuse  ambient 
stray  light  is  tolerable  and  will  not  affect  the  a.c.  drop  signal. 

The  gain  of  the  a.c.  amplifier  should  be  adjusted  here  or  in  the  attenuator 
of  the  interface  receive  electronics  so  that  a  3.6mm  diameter  test  drop  gives 
a  pulse  of  3.6V.  These  drops  can  be  simply  produced  by  a  slow  flow  of  water 
through  a  hyperdermic  needle. 

(c)  Interface  electronics 

The  detailed  circuit  is  shown  in  Figure  6  with  the  block  diagram  and 
timing  sequence  in  Figure  3. 

Considering  the  timing  sequence  we  see  that  when  the  comparator  resets, 
the  peak  detector  will  hold  the  maximum  value  of  the  pulse,  and  this  resetting 
triggers  a  monostable  which  initiates  an  analogue  to  digital  conversion  of 
the  value  held  by  the  peak  detector.  When  the  conversion  is  complete  (about 
20psecs)  the  end  of  conversion  (EOC)  pulse  fires  a  monostable  which  closes  a 
switch  thus  discharging  the  peak  detector  capacitor  and  resetting  it.  The 
eight  bit  output  of  the  A/D  is  connected  to  channel  A  (Port  4  in  the  program) 
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of  the  micro,  whereas  the  end  of  convert  output  of  the  A/D  is  connected  to  the 
most  significant  bit  (MSB)  of  channel  B  (Port  5  in  the  program).  Thus  as  we 
shall  see  in  the  next  section  the  program  sits  in  an  idling  loop  continuously 
monitoring  the  MSB  of  Port  5i  on  detection  of  an  EOC  pulse,  it  reads  Port  4 
which  contains  the  pulse  height  information,  and  then  reads  Ports  8  and  9 
which  contain  the  time  elapsed  since  the  last  pulse  (see  Section  5a). 

The  computer  then  waits  for  another  EOC  pulse. 

Within  30psecs  of  detecting  the  end  of  one  pulse  the  hardware  is  ready 
for  another  pulse.  In  Section  4  it  was  calculated  that  the  minimum  pulse 
separation  even  in  30ms  ^  wind  should  be  140ysecs. 

*  The  531  and  the  JFET  together  form  a  peak  detector  circuit. 

The  lOOOpF  capacitor  after  the  531  is  charged  through  the  diode  for  positive 
going  input  signals  only,  the  JFET  acts  as  a  simple  follower  to  maintain 
high  input  impedance  so  that  once  charged  the  lOOOpF  discharges  only  very 
slowly. 

The  circuit  shown  in  Figure  6  has  two  possible  signal  inputs.  The 
signal  from  the  line  driver  of  the  photodiode  amplifier  can  be  connected 
straight  into  the  peak  detector  for  real  time  analysis;  or  the  signal  from 
the  tape  recorded  raindrop  data  acquired  in  Virginia  can  be  amplified  before 
it  is  routed  to  the  peak  detector.  In  noisy  environments  it  may  be 
necessary  to  use  a  differential  amplifier  to  reject  any  noise  picked  up 
on  both  leads  of  the  twisted  pair  from  the  distrometer,  but  this  has  not  yet 
been  found  to  be  necessary.  The  eight  L.E.D.'s  display  the  latest  value  of 
the  A/D  converter  and  form  a  useful  indication  that  drops  are  being  detected. 
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6  Software  For.  Z80  Eased  Microcomputer 

In  the  present  configuration  of  the  microcomputer,  the  first  8k  of 
memory  are  used  for  'Basic'  with  a  small  amount  of  this  space  reserved  for 
user  written  machine  code  subroutines.  The  memory  from  8k  to  20k,  that  is 
12k  eight  bit  words,  is  used  for  storing  the  data  acquired  by  the  subroutine 
in  real  time,  it  can  be  subsequently  accessed  by  the  basic  program  using  the 
'peek'  statement. 

a)  Z80  Machine  Code  Subroutine 


Qn  receipt  of  a  pulse  from  the  distrometer  we  require  to  store  its 
amplitude,  the  time  since  the  last  pulse  in  quarter  millisecond  units,  and 
occasionally  we  require  the  absolute  time  to  be  recorded  to  enable  the  absolute 
value  of  the  concentrations  of  drops  to  be  calculated.  As  described  previously 

the  memory  from  8k  to  20k  (2000  to  4FFF  in  hex)  is  available  for  this 
data.  This  memory  has  been  split  up  into  256  byte  'pages',  each 
page  starts  with  a  header  comprising  two  zeros  and  then  the  absolute  time  in 
minutes  and  seconds  read  from  the  CTC  (timer)  board  on  the  microcomputer. 

The  data  from  each  drop  then  fills  the  remaining  252  bytes  on  the  page.  Each 
pulse  is  recorded  in  the  format  :  time  in  quarter  milliseconds,  quarter 
millisecond  overflow,  pulse  height.  Thus  on  the  page  we  can  store  252/3  -  84 
pulses,  or  a  maximum  of  42  drops  assuming  there  are  no  drops  giving  single  pulses 
which  will  be  rejected  by  the  basic  pulse  sorting  program.  For  12K  of 
memory,  we  have  48  'pages'  or  a  possible  total  of  48x42*2016  raindrops. 

The  machine  code  subroutine  which  effects  this  data  format  is  shown 
in  Figure  1  The  handshaking  comprises  the  section  of  program  here 


address  OC94  to  OCAD  where  the  program  idles  until  a  transition  from  1  to 
zero  arrives  on  port  5  bit  7  (the  end  of  convert  pulse  from  the  A/D),it 
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then  stores  the  eight  bit  values  of  ports  8,  9  and  4  (16-bits  of  elapsed 
time  in  quarter  milliseconds  and  port  4, the  A/D  value)  in  sequential  memory 
locations  using  the  HL  register  pair  as  a  pointer  for  the  memory  address. 

Turning  now  to  the  other  steps  in  the  program  we  see  that  if  the  HL 
register  pair  hol<k  the  pointer  address  for  storing  the  data  then  as  we 
increment  HL  whenever  H  changes  we  have  just  filled  up  one  page  of  data. 

Register  D  holds  the  current  value  of  HL  and  after  incrementing  HL  we  test 

to  see  if  H^D,  if  this  is  so  a  page  has  been  filled.  Two  courses  of  action 
are  now  possible.  C  is  initialised  as  the  maximum  value  of  H  we  want,  so  if 
H*C  we  have  enough  data  and  we  return  to  the  main  program;  otherwise  we 
initialise  a  new  page  and  wait  for  some  new  pulses. 

Looking  at  the  initialisation  in  detail.  At  the  start  HL  is  set  to  2000 

and  then  in  this  program  C  (*Hmax)  is  set  to  34  for  16^^  pages  of  data. 

C  is  usually  set  by  'poking'  the  value  into  location  OC84  in  the  main  program. 
We  now  enter  the  loop  to  initialise  a  new  page,  and  set  D«H,  and  write  0,0, 
time  in  minutes  (port  11)  and  in  seconds  (port  10),  and  enter  the  pulse 
handshaking  loop.  After  detecting  each  pulse  and  writing  three  bytes 
into  memory  we  check  whether  the  page  is  full.  If  D»H  we  go  and  look  for 
another  pulse,  but  if  D^H  we  check  whether  the  new  H-C  (our  value  of  Hmax) 
in  which  case  we  return  to  the  main  program.  If  H?*C  then  we  go  and 
initialise  a  new  page  at  OC83.  However  to  indicate  to  the  operator  that 
data  is  being  gathered  and  reassure  him  that  the  computer  really  is  filling 
its  memory  with  data^before  going  to  OC83  we  jump  to  a  subroutine  in  the 
computer  monitor  (instruction  DF  66)  which  simply  transfers  the  contents  of 
the  HL  register  to  the  TV  monitor.  As  each  page  is  filled  up,  the  address 
of  each  new  page  of  data  is  displayed  for  the  operator  to  see. 


b  Basic  Program 


The  basic  program  developed  for  sorting  the  pulses  acquired  in  the  machine 
code  subroutine  into  spectra  and  outputting  the  data  on  to  a  small  printer  is 
shown  in  Figure  3.  It  should  be  stressed  that  this  program  is  not  optimised 

and  provides  many  output  statements  for  monitoring  the  performance  of  the 
distrometer  which  would  not  be  required  for  a  non-prototype  instrument.-  The 
program  outputs  spectra  calculated  by  the  two  independent  methods  so  as  to  enable 
a  comparison  to  be  made.  These  methods  are  the  'flux'  method  with  an  assumed 
drop  terminal  velocity  and  the  'time  of  flight'  method  whereby  no  velocity 
information  is  required  but  the  spectra  are  derived  from  knowing  the  fraction  of 
the  total  time  that  a  drop  of  a  certain  size  is  within  the  sample  volume. 

The  program  is  fairly  straightforward: 

1  15-70  set  up  the  ports  and  call  the  m/c  code  subroutine 
l  90-150  zero  various  arrays  used  for  storing  and  counting  drops 
1  155-175  output  the  amplitudes  and  times  of  the  pulses  on  the  first 
page  to  the  printer 

l  205-  For  each  new  page  of  data  calculate  the  amplitudes  and  elapsed 
time  between  pulses,  accept  pulse  pairs  with  time  separation 
.  between  1  and  20msecs  (£330-332)  and  amplitudes  agreeing  to 

within  20%  (1  338).  Increment  the  appropriate  bin  counter  (1367) 

1  500  output  the  spectra  calculated  by  the  time  of  flight  and  flux  method 
for  this  one  page. 

1  655  calculate  the  rainfall  rate. 

1  678  Get  a  new  page  of  data  (at  line  205). 

1  680  When  all  the  pages  of  data  analysed^print  out  the  cumulative  spectra 
and  the  equivalent  Marshall  Palmer  spectrum. 


The  disdrometer  is  designed  to  be  simple  to  use  and  not  to  require  any 
elaborate  alignment  procedure.  Providing  the  frame  on  which  the  collimating 
and  receiving  optics  are  mounted  is  reasonably  rigid  the  instrument  should 
only  need  aligning  once.  The  photodiode  may  be  moved  in  a  vertical  plane 
as  shown  in  the  engineering  drawings,  and  its  position  should  be  adjusted 
until  the  light  intensity  as  measured  by  the  DC  monitor  level  output 
(Section  5b)  is  a  maximum.  This  value  should  be  about  6V.  At  this  stage 
all  the  light  of  the  annulus  should  be  focussed  on  to  the  photodiode. 

A  further  simple  check  is  to  use  drops  produced  by  a  slow  flow  of  water 
through  a  hyperdermic  needle  to  check  that  the  pulse  pairs  are  of  equal  size 
and  do  not  vary  with  position  around  the  annulus,  or  along  the  beam.  As 
mentioned  in  Section  5b,  the  gain  of  the  preamplifier  (or  more  conveniently 
the  attentuation  between  the  preamplifier  and  the  interface  electronics) 
should  be  adjusted  so  that  a  3.6  mm  drop  produces  a  3.6V  pulse.  This  can 
be  verified  either  using  the  LED  display  of  the  interface  electronics  or  by 
running  the  program  and  checking  that  all  the  pulses  are  sorted  in  the  cor¬ 
rect  bin.  Once  set  up  the  disdrometer  should  not  go  out  of  alignment  but 
the  above  test  using  3.6  mm  drops  should  always  be  used  as  a  check. 

b)  Results 

Figures  9  and  10  display  some  spectra  obtained  in  Virginia,  USA  and 
Manchester,  UK.  We  can  see  in  Figure  9  that  there  is  reasonable  agreement 
between  the  spectra  calculated  for  the  88  second  period  by  the  conventional 
droplet  flux  method  and  the  novel  time  of  flight  criterion  unique  to  this 
disdrometer.  The  counts  displayed  in  bin  1  (200-400um  radius)  may  be  unreliable. 
Both  spectra  show  less  small  drops  than  the  equivalent  Marshall  Palmer.  To 
calculate  the  significance  of  these  spectra  the  second  column  labelled 
"cumulative  number"  in  the  Table  gives  the  total  number  of  drops  (n)  recorded 
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in  each  bin.  From  Poisson  statistics  the  fractional  error  in  the  concentration 
for  that  bin  is  l//n.  Consequently,  we  can  see  that  although  the  spectra 
calculated  by  the  two  methods  agree  for  the  larger  drops,  the  flux  method 
is  recording  more  drops  in  bins  2,  3  and  4.  This  is  to  be  expected  as  there 
was  a  light  wind  of  a  few  metres  a  second  which  would  increase  the  apparent 
flux  of  smaller  drops  more  than  the  larger  ones.  The  time  of  flight  method 
is  of  course  unaffected  by  any  wind. 

The  spectra  shown  in  figure  10  was  taken  in  very  light  drizzle.  Again 
the  wind  speed  was  a  few  metres  a  second  and  we  see  that  far  more  small  drop¬ 
lets  are  recorded  using  the  flux  method  than  the  time  of  flight  principle. 
However,  even  the  time  of  flight  method  is  recording  more  drops  than  predicted 
by  the  Marshall  Palmer  for  the  equivalent  rainfall  rate.  Although  the  statistics 
in  bin  3  (total  drop  count  9,  error  33%)  are  not  good,  it  may  well  be  that 
the  often  quoted  statement  that  the  Marshall  Palmer  overestimates  the  number 
of  small  droplets  may  not  be  true  for  the  fine  light  drizzle  so  prevalent 
in  the  United  Kingdom. 
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8  Future  Work 


The  final  Figure  11  shows  a  photograph  of  the  disdrometer  on  a  dexion 
framework  for  field  operation.  Versions  developed  during  the  past  twelve 
months  were  rather  larger  using  collimating  lenses  of  focal  length  1m  and 
then  5ccm,  rather  than  the  25cm  lens  used  in  the  version  in  this  photograph 
and  in  the  engineering  drawings.  It  would  be  difficult  to  reduce  the  size 
further.  As  has  been  described  in  this  report  the  optical  set-up,  the  elec¬ 
tronics  and  the  software  have  been  developed  and  are  working  satisfactorily. 
Results  show  that  the  new  method  of  calculating  raindrop  concentrations  using 
the  time  of  flight  of  the  drops  across  the  sample  volume  works  well.  In  windy 
conditions  the  conventional  flux  method  overestimates  the  concentration  of  drops. 
It  is  worth  noting  that  a  method  which  relies  on  measuring  the  flux  across  a 
horizontal  surface,  such  as  the  Knollenberg  prober,,  will  have  problems  due  to 
the  gusty  nature  of  the  wind.  The  small  drops  in  windy  conditions  do  not  fall 
at  a  constant  angle  but  are,  on  occasion  even  observed  rising  upwards. 

More  work  and  specifically  more  experience  of  using  the  present  disdrometer 
in  many  different  types  of  naturally  occuring  rainfall  are  needed  before  the 
instrument  can  be  deemed  truly  operational. 

Some  of  the  remaining  tasks  are  listed  below. 

a)  Field  trials  comparing  the  rainfall  rate  with  a  tipping  bucket  raingauge 
and  if  possible  another  disdrometer. 

b)  Operation  in  high  wind  speeds.  A  comparison  of  the  two  methods  of  calcu¬ 
lating  the  raindrop  spectra  whilst  measuring  the  windspeed.  An  examination 
of  the  statistics  of  the  transit  times  of  the  drops  within  a  certain  bin 
size  to  see  if  they  agree  with  theory,  and  are  modified  by  the  ambient  wind 
as  predicted. 
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c)  Drops  landing  on  the  lenses.  The  problem  of  drops  landing  on  the  lenses 
in  the  vicinity  of  the  annular  mask  is  not  as  severe  as  with  scattering 
instruments.  The  direct  transmission  of  light  through  a  drop  on  the  lens 
is  not  greatly  attenuated.  However,  in  Virginia  a  rather  primitive  heater 
and  blower  were  installed  to  keep  the  lenses  warm.  This  arrangement  worked 
well  but  the  wind  speeds  were  never  excessive. 

If  we  observe  that  the  light  to  be  collected  by  the  lenses  is  only  in 
a  ring  100pm  wide  but  that  the  raindrops  hitting  the  lens  will  have  had 
trajectories  at  various  angles,  it  is  obvious  that  a  serrated  shield  of 
multiple  discs  rather  like  a  high  voltage  insulator  mounted  on  the  centre 
of  the  lens  will  allow  the  parallel  annular  beam  of  light  to  fall  on  the 
lens,  but  make  it  virtually  impossible  for  the  raindrops  and  any  splash 
products  to  hit  the  lenses.  Also  any  arrangements  for  blowing  air  onto 
the  lenses  need  be  less  powerful  because  the  lenses  themselves  will  be 
less  well  ventilated  by  the  ambient  saturated  air. 

d)  Splashing.  In  heavy  rain  the  disdrometer  will  also  measure  splash  products 
passing  through  the  beam.  Both  ground-based  commerical  instruments  suffer 
from  this  problem.  However,  the  device  (as  shown  in  the  photograph)  can 

be  mounted  some  distance  above  the  ground  to  stop  drops  bouncing  from  the 
ground.  Sponge  shields  are  placed  on  the  upper  surface  of  the  transmitting 
and  receiving  optics  (not  shown  in  the  figure)  and  in  the  light  rains  so 
far  encountered  they  have  worked  well.  One  can  envisage  two  disdrometers 
at  right  angles  but  with  intersecting  sample  volumes;  by  accepting  only 
drops  observed  by  both  instruments,  the  drops  in  this  intersecting  volume 
are  measured,  and  are  supposedly  far  from  the  solid  supports  and  therefore 
with  no  splash  products.  However,  an  examination  of  the  statistics  and 
especially  the  pulse  arrival  rate  Bhows  that  this  is  not  a  very  practical 
method. 

e)  Two  sample  volumes.  The  present  device  with  a  sample  volume  17cm  long 

and  3cm  acroas  is  a  compromise  between  the  diffraction  limit  for  small  drops 
and  the  need  to  measure  a  significant  number  of  the  rarer  large  drops. 
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However,  to  attain  a  meaningful  sample  of  all  sizes  of  drops  rapidly,  two 
versions  could  be  run  together,  a  small  sample  volume  for  the  small  drops  and 
a  bigger  one  (up  to  lm  50cm  long)  for  the  drops  larger  than  1mm  diameter. 

f)  The  present  disdrometer  is  portable  but  requires  considerable  ancillary 
equipment  consisting  of  the  computer  and  keyboard,  TV  monitor  and  printer. 

A  more  ideal  solution  would  be  to  burn  the  debugged  and  finalised  program 
on  to  an  EPROM  and  output  the  data  on  to  a  digital  mini-cassette.  The  elec¬ 
tronics  would  then  all  fit  in  the  present  box  containing  the  computer  interface. 

g)  Lower  power  version.  Most  of  the  power  consumed  by  the  device  is  in  the 
illumination  (12V).  For  remote  sites  a  low  power  version  could  be  developed. 
Firstly,  less  illumination  could  be  used  but  more  amplification  in  the 
electronics.  Secondly,  a  tipping  bucket  rain-gauge  could  be  used  to  switch 
the  light  on  during  rain. 
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TABLE  1 


Minimum  and  Maximum  Pulse  Separation  Times 


DROP  DIAMETER 

MIN.  TRANSIT 
DISTANCE 
(mm) 

TERMINAL  VELOCITY  I 

30  ms- 

1  WIND 

V  1 
ms 

MIN 

msecs 

MAX 

msecs 

MIN 

msecs 

MAX 

msecs 

300pm 

4.2 

1.15 

3.6 

26 

0.14 

i 

500pm 

5.5 

2 

2.7 

15 

0.18 

i 

1  mm 

7.5 

4 

1.9 

7.5 

0.25 

i 

2  mm 

11 

6.5 

1.7 

4. 6 

0.36 

i 

3  mm 

13 

8.0 

1.6 

3.8 

0.43 

i 

4  mm 

15 

8.8 

1.7 

3.4 

0.5 

i 

The  maximum  transit  distance  is  always  30  mm 


TABLE  2 


Probability  of  one  drop  (f(l))  and  two  drops  (f  (2))  in  the  sample 
-3  3 

volume  (0.12  10  m  )  assuming  a  Marshall  Palmer  raindrop  distribution 


RAINFALL  RATE 

,  -1 
mm  nr 

DROPS 

>  500um 

DROPS 

>  300pm 

n  (m  3) 

f  (1) 

f  (2) 

n  (m  3) 

f(l) 

f  (2) 

1 

230 

.02 

.0002 

540 

.06 

.002 

3 

450 

.05 

.001 

900 

.1 

.005 

10 

850 

.1 

.005 

1400 

.17 

.015 

FIGURii  3 


b) Timing  Diagraa 


IDlif  iar 


Car.^irator 


I'lonoatable 


"i-Out* c 


A/D  ROC 


Monostflble  ? 


FISUHiS 


External  electronics  for  the  timer  counter  board. 


^  l/rrw  Q-c-txrcC  -AjCyU/yfJZc^ 

s)  c^vO  —  I Nf\rr  4  kW>  ^  b-rr^k  ucc*v5^  'ZS'S 
(Potcr  i)  C-M  —  \U<KT  e*0  <£rrx  utirW"  2$S 


(^PORT  io)  cPvH 
(^PORT  llj  c&3 


I  NPcr  1  Ri 

input  cS.,3 


^£rr»a  vXXT'Vjf  ^  0 
^ur*a  vcevaf  ^  O 


O'-^aet’  i5S  A 

-tArl-Wj  iiS  X.  -i.  , 

rrtajJ-  Aj//vcrr 

/rea»(t  *e\r&n£  /vc^^v  - 


4-001 


N  tj  -r  V 


r  1* 

*  — 

^  2  X 

CrU-suA  4-IH  =  2  Kv 
■J  nki . 


Sock 6 t  c 


Function  P‘>U 


0  TAl  ^  10  <—  V074k»  _in 

O  °/p  T  -v 

“  )  LlVK 

4.  Tfll^  1 1  * 


4094  k*  lkl 

L  L 

h-^5V% 

?3JK- 


SVsuts  I  k> 


FUMCmv  P«V 

%  Ta.<;  n 

*  %  5  ' 

3  mu;  i3k 


PCBT  ^rTTlM^  ~  \jsdj.  ^><frht  -un  JCCW^ltv  STMtlg . 

cut  r(  m  f  our  a ss  •  out  9;  m  :  out  9,a*.r 

OuT  ic  m7  :  Out  ic,  40  Out  It.  U7  OUT  U.  <0 


FIGlhtE  5  PRE-AMi  LIFIEH  AND  LIME  DRIVER 


$ 

3  *  3 

4  ^ 

I  5  |g 

o  i  a 
NR 


AOORESS  CONTEXTS 


0C80  21  ] 

1 

*\ 

0C81  00 

i.0  m.  9  iooo 

0C82  20  J 

1 

0C83  0E  1 

l  LO  C  34 

0C94  34 

0C85  54 

LO  0,  W 

/ 

0C36  3E  1 

1  LO  A,  w  0  ^ 

0C87  00  J 

0C88  77 

(ml)  ♦—  A 

0C89  23 

INC  hl 

0C8A  77 

(hl)  A 

0C8B  23 

INC  hl 

0C8C  DB  1 

[  hv  A,  08 

0C8D  0B  ' 

t  * 

0C8E  77 

*  A 

0C8F  23 

INC  Hu 

0C90  DB  1 

[  IN  A  OA 

0C91  0A  ! 

l 

0C92  77 

(}tL)«-  A 

0C93  23 

INC  HL  J 

0C94  DB 

\  Itv  A. 5 

0C95  05 

J 

0C96  CB  , 

1  Bit  *7,  A 

l 

0C97  7F 

1 

j 

) 


0C98  CA 
0C99  94 
0C9A  0C 
0C9B  DB 
0C9C  05 
0C9D  CB 
0C9E  7 F 
0C9F  C2 
0CA0  9B 
0CA1  0C 
0CA2  DB  I 
0CA3  08  J 
0CA4  77 
0CA5  23 
0CA6  DB 
0CA7  09 
0CA8  77 
0CA9  23 
0CAA  DB 
0CAB  04 
0CAC  77 
0CAB  23 
0CAE  7A 
0CAF  BC 
0CB0  CA 
0CB1  94 
0CB2  ec 
0CB3  79 
0CBA  BC 
0CB5  C8 
0CB6  DF 
0CB7  66 
0CB8  C3 
0CB9  93 
0CBA  BC 
0CBB  C9 
0CBC  C9 
0CBD  C9 


|  TP  i  oc<t* 
1 


IN  A ,S 
Bit  1,  A 

TP  Ml  OOB 


IN  A,  8 


RtTUAN 


FIGURE  Tt:  280  MACHINE  CODE 
SUBROUTINE 

Initialise  constants 
(C  poked  down  from  main  program) 

Initialise  'page'  of  memory  with  0,0, 

time  in  minutes  (ch11=0B),  and  seconds  (ch10=0A) 


Wait  here  until  MSB  port  5  *  1 


Wait  here  until  MSB  port  5=0 
i.e.  exit  on  negative  edge  on  MSB 


Store  channels  3,9,  and  4  in 
sequential  memory  locations. 
The  HL  pointer  is  incremented 
after  each  port  is  read. 


If  D  «  H,  page  is  not  full 
go  and  wait  for  more  pulses 

If  G  (“H  v)  =  H,  enough  data 
return  to  the  main  program 

Page  full.  Output  (HL)  on  to  monitor, 
then  go  and  initialise  and  fill  another 
page 


BASIC  PROGRAM 


FIGURE  a  a 


] 


3 


1 5  OUT8 » 1 1 7 :0UT8 r  255 : OUT9  r 1 17 : OUT9  , 255 
28  OUT10> 117*OUT10»60:OUT11, 117:OUTil ,60 

25  0UT6 1 255  s  0UT6 » 255 : 0UT7  r 255 : OUT7 » 255 

26  INPUT"NUMBER  OF  PAGES" »Y: Z=Y+64 

27  INPUT"NUMBER  OF  BINS" ?N:PRINT 
50  IF<  Z>80 >G0TO26 
70  POKE3204»ZsDOKE4100»3200sA=USR< 0 )  J  Call 
80  B=INP< 11 >SC=INP< 10  >:A=16384+< 256*< Z-64 )  ) 
90  POKE  At 0 SPOKE  A+1,0:POKE  A+2,B:POKE  A+3,C 

96  DIM  DD(N)»FFCN)fD(N)fF(N),R<N)fS<N) 

97  DIM  RAD<  N),GF<  N  >»VEL<  N) 

98  DIM  Xl(N),X2(N>fNUM(N)rE(N>fRl(N>fR2(N) 
100  F0RK»1T0N:£<  K  >=0sDD<  K  >=0:FF<  K  >=0:NEXTK 


initialise  Ports 


SET  Of  HA*  HO  PA^tf 


SutMvnut . 

\  TBRmwMt  Tme  on 


2  ERO  ARRAYS 


RETURN 


150  WIDTH90  s P0KE4 1 63  ,  (  INT<  90/14  )-l  )*14 
155  I*1:A*16388 

160  B*PEEK<  A)*C=PEEK<  A+l  )*  W=128-PEEK<  A+2  ) 
165  A=A+3* 1*1+1 : IF< I >84  )THEN  PRINT JGOTO200 
170  T=0.25*( ( BB-B  )+256*< CC-C ) ) 

175  BB*BsCC*CsPRINTT;W»sGOTO160 


CAlcutate  And  output  volts 

Amo  elappeo  tib*  r>r  f<rst  PA$f 
.  of  data. 


200  AVR*0 

205  FORM* 1TOY *SUM=0:SSU*0 

222  FORK*lTON s RAD<  K  )*(  (  100*N  >+50  >*0.001 

224  GF<  K  >*<  C  <  RAD<  K  )t3  >#1 . 33*3 .  14 16  >/2 . 9  >t0  . 382 

226  VELC  K  )*9.32032*<  l-<  l/(  EXP<  GFC  K  ) )  )  ) ) 

310  D<K)*0rF(  K  )*0:NEXTKSI*1 

314  A*1 6388+256*< M-l  >*OM=PEEK(  A- 2  ):OS=PEEK<  A-l  ) 
317  BB*PEEKC  A  )sCC=PEEN<  A+l  ):  W=128-PEEK<  A+2  ) 

320  A*A+3 *  1*1+1 »IF< I >84  )GOTO410 
323  B*PEEK(  A )*C=PEEN<  A+l  )*V=128-PEEK<  A+2  >*A=A+3 
326  1*1+1 s IF(  I >84  >GOTO410 
328  T*0. 25*<  <  BB-B  )+256*<  CC-C  ) ) 

330  IF< T<1  )GOTO400 
332  IF<  T>20  )GOTO400 
334  VW*V+W:R*V:IF<  UU>V  )THENR=W 
330  IF(  ABS(  V- W  >  > . ?*R  )GOTO350 
340  K*INT< 0.005*VW*N > 

342  E(K)*E<  K )+ 1  s D<  K  )*D<  K  )+l:F<K  )=F(  K  )+TsG0T0317 
350  NB*PEEK( A  )sNC=PEEKC  A+l >sNV=128-PEEN<  A+2  ) 

352  A*A+3* 1*1+1: IF< I >84  )GOTO410 
355  T*0 . 25* < ( BB-NB  )+256*<  CC-NC ) ) 

360  IF<  T< 1  )GOTO380 

362  IF( T>20  )GOTO380 

363  R*NV  s IF<  VV  >NV  )THENR= W 

364  IF( ABS<  NV— W  )> . 9*R  )GOTO380 

365  VW*W+NViK*INT<  0.005*VW*N  ) 

367  E<K  )*E<K)+1*D(K  >*D(  K)+UF<  K  )*F<  N)+T:GOTO400 
380  BB*B *  CC«C  s  VV= V  * B=NB : C*NC s  V=NV : G0TO328 
400  BB*B  *  CC*C  *  UV*U  s  G0T0323 
410  NM*PEEK<  A+2  )*NS*PEEK<  A+3  > 

420  IF<  OM<  NM )THENNM=NM-60 
425  C0T0432 

430  PRINT*  PRINT"PAGE"  M 
432  TT*( OM-NM  )*60+<  OS-NS ) 

435  PRINT"FROM"QM"MINS" IOS"SECS  "? 

440  PRINT"TO"NM"MINS" ?NS"SECS. " t 
445  PRINT"TIME  LAPSE* "TT"8ECS" SPRINT 
448  IF<  M<Y  )COTO480 

450  PRINT  "  RADIUS  No.  CUM. No."? 

460  PRINT"  TIME(m*>  C0NC'N./*t3  "I 
470  PRINT"  CUM. SPECTRUM" 

473  PRINT  SPC<  38  >"  FLUX  T.OF.F"? 

475  PRINT  SFC<7)"FLUY  T.n^.F- 


FOR  EACH  PACE  (h) 


] 


CALCULATE  TERMINAL 

velocities 


FDR  each  Pulcb  Pair. 

CALCULATE,  AMPLITUDE 


> 


J 


Ano  Elapsed  time. 

ACCEPT  TIMES  V  -  '20  nvr«t 

And  pairs  A^rEEin^  ^ 

WITHIN  2.0% 

K  16  Bint  nuhoer. 


INCREMENT  BIN  COUNTS 


LOOK  FOR 
OVERLAPPING  POL.SE 


Pairs 


WRITE  TITLE  FOR 
SPECTRUM 


FIGURE  6 


BASIC  PROGRAM 


WIDTH90:POKE4163 1  <  INK  90/14  >-l  )*14 

F0RK  =  1T0N:R<K  >*D<  K  )*200*17/<  23*TT*VEL<  K>> 

R(K  >*<  INK  R<  K  >*100  >  )/100 

S<  K  )=F<  K  )*17/(  23*0. 12*TT  > 

S<  K  >*<  INK  S<  K  )*100  )  )/100 

DD<  N  >=DD<  K  )+R(  K  )  :FF(  N  )=FF(  K  )+S<  K  ) 

Rl<  K  )=INT<  2000*K/N  )sR2(  K  )*INT<  2000*<  K>1  )/N  ) 
Xl<  K  )-R<  K  )*VEL<  K  )*RAQ<  K  >+3*36/lE4 
SUM*SUM+X1< K) 

X2<  K  >*S<  K  )*<  <  2*RAD<  K  )  )t3  )*VEL<  K  ) 

SSU=SSU+X2<  K )JNEXTN 
IF<M<Y  >G0T0655 
FORK*  1  TON  s  PRI NTK 5  TAB(  4  >  5 
P*=STR*<  R2<  K  )  > 

P**RIGHT*<  P*»<  LEN(  P*  )-l  )  ) 

J**STRS<  Rl<  K  )  )+"-"+P*+"um" 

PRINTJ*5TAB(17)5 
PR INTD<  K  >  5  TAB<  23  )  5  E<  K  )  5  TAB<  31)5 
547  PRINTFt  K)5TAB(41  )5 
PRINTRt  K)5TAB<  51  )5S<K  )5TAB<  61  )5 
PRINK  INK  DD<  K  )/M*100  )  )/1005  TAB<  71  >5 
PRINK  INK  FF<  K  >/M#100  )  )/100  s  NEXTK  SPRINT 
PRINT"SPECTRUM  FOR  PAGE  8  "M 
F0RK=1T0NSPRINTK5TAB(  4 >5 
PRINK  INK  R<  K  )  )  )5TAB<  9  )5 
IFR<  K )*0THEN  PRINT SGOTO600 
FORL*1TO(  10*LOG(  R<  K  )  )  ) 

PRINTCHR*< 42  )5 SNEXTL8PRINT 
NEXTK 

PR=<  INK  <  SUM*1  .33*3.14159  )*100  )  )/100 
RT=<  INK  <  SSU*<  6*3. 14159  )*lE-4  )*100  )  >/100 
AVR*AVR+PR s PRINT :G0T0678 

PRINT-R/RATE  <  tav  FLUX  Method  ):  "PR 5  "MM/hr" 
PRINT"R/RATE  < bv  TIME  OF  FLIGHT  Method  >s  "5 
PRINTRT  5 "MM/hr" s  PRINT 
i  NEXTM 

PRINT "CUMULATIVE  SPECTRUM< I )" 

FORK* 1 TON s PRI NTK 5 TAB< 4 )5 
PRINK  INK  DD<  K  >/(  M-l  )  )  )5TAB<  9  )5 
IFDD<  K )/( M-l  )*0THEN  PRINT8GOTO690 
FORP*lTO<  10#LOG(  DD<  K  )/<  M-l  >  )  > 

I  PRINTCHR*< 42)5 sNEXTPsPRINT 
NEXTK 8 PRINT 

PRINT"CUMULATIVE  SPECTRUM<  ID- 
FORK*  1TONSPRINTK 5  TAB<  4  )5 
PRINK  INK  FF<  K  )/<  M-l  )))5TAB(9)5 
IFFF<  K  )/(  M-l  )«0THEN  PRINT 8G0T069? 

’  FORP*lTO< 10*LOG<  FF( K )/(  M-l ) ) ) 

PRINTCHR«  42  )5  8NEXTPJPRINT 
i  |^XTK 

N0*8000  s  RR* AVR/Y  8  B I N=0 . 28  CN* -0 . 2 1 » PRINT 
PRINT  "AVERAGE  R/RATEs»lNT<RR*100  )/1005 
;  PRINT"MM/hr" 

PRINT  8Pl*RRtCN8P2=Pl*4. 1 
FORK-1TON  8  P3*P2*<  RAW  K  >*2") 

NUM<  K  )«INT<  <  N0*EXPt  -P3  >  >*BIN  >*NEXTK 
I  PRINT"EQUIVALENT  M-P  SPECTRUM" 
F0RK=1T0N8PRINTK!TAB<  4  )5 
i  PRINTNUM<  K)!TAB<9>» 

IFWJM< K  )»0THENPRINT 1GOTO790 

F0RP»1T0< 10*LOG< NUM< K )  )  )8PRINTCHR*<  42  )5 

NEXTPlPRINT 

NEXTK 


FbR  Each  pa<;e  of  data 


CALCULATE  SPECTRUM 


usiwj 


FLUX  A  NO  T>  HE 


op  flight  newop. 


Output  spectrum 


OP  PA^E  OF 


DATA. 


CALCULATE  RAINFALL 
RATES 


NE*T  PA^g . 


OUTPUT  CUMULATIVE 

Spectra  found 

USlN£  FLUX  AND  TIME 
OP  FLIGHT  METHOO . 


Output  Equivalent 

MARSHALL  -  PALMER 


SPECTRUM. 


FIGURE  9  RAl.iDROF  SPECTRUM  IN  ZIHGII.IA.  TOTAL  3ARPL3  88  dESUNDS 


RADIUS 

No . 

CUM .  No . 

TIME'  ms  ) 

CONC'N. 
FLUX  T 

/m*3 
•  OF.F 

CUM. SPECTRUM 
FLUX  T.OF.F 

1 

200-400UM 

2 

26 

17.75 

43.14 

18.22 

38.42 

14.67 

2 

400-680um 

5 

146 

48.25 

63.16 

49.53 

135.41 

87.11 

3 

600-800UM 

12 

136 

71.75 

108.65 

73.65 

97.05 

61.81 

4 

800-1000um 

6 

99 

22.75 

42.98 

23.35 

51.21 

38.17 

5 

1000-1200UM 

5 

60 

26.5 

30.05 

27.2 

25.46 

21.58 

6 

1200-1400u« 

2 

26 

8.25 

10.48 

8.4’6 

9.72 

7.22 

7 

1400- 1600un 

0 

6 

0 

0 

0 

1.77 

1.5 

8 

1600-l800un 

0 

o 

0 

0 

0 

.85 

.84 

9 

1800-2000UM 

1 

r> 

4.25 

3.99 

4.36 

.44 

.59 

10 

2000-2200u« 

0 

0 

0 

0 

0 

0 

0 

SPECTRUM 

1  43 

2  63 

3  108 

4  42 

5  30 

6  10 

7  0 

8  0 

9  3 

10  0 


FOR  PAGE*  15  {last  6  seconds  of  sample  period) 

%X%%%\ t*********************************** 
********************************************** 
************************************* 
********************************** 

*********************** 


************* 


(units  are  concentratior 
in  number  per  cubic 
metre) 


CUMULATIVE  SPECTRUM< I  ) 

1  38  ************************************ 

2  135  ************************************************* 

3  87  ******************************************** 

4  51  *************************************** 

5  25  ******************************** 

6  9  ********************** 

7  1  ***** 

a  0  * 

9  0  * 

18  0 


Flux  method  : 
Spectrum  for 
88  second 
sample 


CUMULATIVE  SPECTRUM<  II  ) 

1  14  ************************** 

2  87  ******************************************** 

3  61  ***************************************** 

4  38  ************************************ 

5  21  ****************************** 

6  7  ******************* 

7  1  **** 

8  8  * 

9  0  * 

10  0 

AVERAGE  R/RATEi  1.04  Mfi/hr 

EQUIVALENT  M-F'  SPECTRUM 

1  472  ************************************************************* 

2  209  ***************************************************** 

3  93  ********************************************* 

4  41  ************************************* 

5  18  **************************** 

6  8  ******************** 

7  3  ********** 

8  1  * 

9  0 

10  0 
Ok 


Time  of  flight 
method,  spectrum 
for  same  sample 
period 


i 

-t 


A 


FIGUR3  10  rtAINJROF  5P5CTRUM  IN  MAiiCHESTiR  TOTAL  SAMPLE  .45  SECONDS 


T". 


RADIUS 

No. 

CUM. No. 

TIME( ms > 

CONC'N. 

/m+ 3 

CUM. SPECTRUM 

FLUX  T 

.OF.F 

FLUX 

T . OF.F 

1 

200-400UM 

7 

118 

75 

226.49 

115.48 

343.78 

151. 

n 

400-600UM 

12 

105 

91.75 

227.4 

141.28 

167.28 

103. 

3 

600-B00UM 

0 

9 

0 

0 

0 

10.22 

6.62 

4 

800-1000UM 

0 

3 

0 

0 

0 

3.58 

4.23 

5 

1000-1200UM 

0 

0 

0 

0 

0 

0 

0 

6 

1200-1400um 

0 

0 

0 

0 

0 

0 

0 

7 

1400-1600um 

0 

0 

0 

0 

0 

0 

0 

8 

1600-1800UM 

0 

0 

0 

0 

0 

0 

0 

9 

1800-2000UM 

0 

0 

0 

0 

0 

0 

0 

10 

2000-2200um 

0 

0 

0 

0 

0 

0 

0 

SPECTRUM  FOR  PAGE*  12 

1 

226  ****************************************************** 

2 

227  ****************************************************** 

3 

0 

4  ®  (last  4  seconds  of  sample  period) 

5  0 

6  0 

7  0 

8  0 

9  0 

10  0 

R/RATE  (bv  FLUX  Method )s  .11  MM/hr 

R/RATE  <  bv  TIME  OF  FLIGHT  Method):  .07  MM/hr 


CUMULATIVE  SPECTRUM 

1  343  ********************************************************** 

2  1 67  *************************************************** 

3  10  *********************** 

4  3  ************ 

5  0 

6  0 

7  g  Flux  method;  spectrum  for  45  second  period. 

®  0  (time  of-  flight  spectrum  is  given  in  the  last  column 

9  0 

10  g  of  the  table  above) 


AVERAGE  R/RATE*  .12  MM/hr 


EQUIVALENT  M-P  SPECTRUM 

1  239  ****************************************************** 

2  67  ****************************************** 

3  18  **************************** 

4  5  **************** 

5  1  * 

6  0 

7  0 

e  e 


END 

DATE 

FILMED 

9-8 


DT  C 


